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We review the result of SUSY parameter fits based on frequentist analyses of experimental constraints from 
electroweak precision data, (g — 2) M , B physics and cosmological data. We investigate the parameters of the 
constrained MSSM (CMSSM) with universal soft supersymmetry-breaking mass parameters, and a model with 
common non-universal Higgs mass parameters in the superpotential (NUHM1). Shown are the results for the 
SUSY and Higgs spectrum of the models. Many sparticle masses are highly correlated in both the CMSSM and 
NUHM1, and parts of the regions preferred at the 68% C.L. are accessible to early LHC running. The best-fit 
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1. INTRODUCTION 

Supersymmetry (SUSY) is one of the favored 
ideas for physics beyond the Standard Model 
(SM) that may soon be explored at the Large 
Hadron Collider (LHC). Here we review the re- 
sults from frequentist analyses [3(2] of the pa- 
rameter spaces of the constrained minimal su- 
persymmetric extension of the Standard Model 
(CMSSM) — in which the soft SUSY-breaking 
scalar and gaugino masses are each constrained to 
universal values mo and 77^ / 2 , respectively — and 
the NUHM1 — in which the soft SUSY-breaking 
contributions to the Higgs masses are allowed to 
have a different but common value. Both mod- 
els have a common trilinear coupling Aq at the 
GUT scale and tan/3 (the ratio of the two vac- 
uum expectation values) as a low-energy input. 
An overview about the existing literature on this 
subject can be found in Ref. pQ. 

The fit includes precision electroweak data, 
the anomalous magnetic moment of the muon, 
(g — 2) M , i?-physics observables (such as rates for 
BR(6 -> sj) and BK(B U —} tv t ), and the upper 
limit on BK(B S —> /i + /j - )), iiT-physics observ- 
ables, the bound on the lightest MSSM Higgs bo- 
son mass, Mh, and the cold dark matter (CDM) 
density inferred from astrophysical and cosmolog- 



ical datqj, assuming that this is dominated by the 
relic density of the lightest neutralino, fl x h 2 . 

2. THE MASTERCODE 

We define a global x 2 likelihood function, which 
combines all theoretical predictions with experi- 
mental constraints: 
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Here N is the number of observables studied, 
Ci represents an experimentally measured value 
(constraint), and each Pi defines a prediction 
for the corresponding constraint that depends 
on the supersymmetric parameters. The exper- 
imental uncertainty, ct(Cj), of each measurement 
is taken to be both statistically and systemati- 
cally independent of the corresponding theoreti- 



x Wc did not include the constraint imposed by the exper- 
imental upper limit on the spin-independent DM scatter- 
ing cross section crjj 1 , which is subject to astrophysical and 
hadronic uncertainties. 



cal uncertainty, a {Pi), in its prediction. We de- 
note by x 2 {Mh) (see below) and x 2 (BR(-B s -> 
/i/i)) the x 2 contributions from two measure- 
ments for which only one-sided bounds are avail- 
able so far. Similarly, we include the lower 
limits from the direct searches for SUSY parti- 
cles at LEP [3] as one-sided limits, denoted by 
"X 2 (SUSY search limits)" in Eq. Q. 

In order to derive a reliable theoretical predic- 
tion of an observable (denoted as Pi in Eq. ([T])) 
and to obtain a correspondingly small theory (in- 
trinsic) uncertainty (denoted as <r(P;) in Eq. ([1])) 
it is crucial to have higher-order corrections (be- 
yond the one- loop level) at hand. As an ex- 
ample, including SM corrections up to four-loop 
and MSSM corrections up to the two-loop level, 
the theory uncertainties of the W boson mass, 
Myy, and the effective leptonic weak mixing an- 
gle, sin (9 ff are now at the level of (depending in 
the MSSM on the SUSY mass scale) 
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Other examples are BR(6 — > 57), where in the 
SM the intrinsic uncertainty is now at the level of 
~ 0.27 x 10 -4 [8] (adding errors in quadrature) , 
or {g — 2) M with a theory uncertainty of less than 
4.9 x 10 -10 [5]. Many calculations and tools used 
in our analyses were presented at "Loops & Legs" 
conferences over the past decade, see Ref. |10j . 

The main computer code for our evaluations 
is the MasterCode [T] 12] ITlTfl"3] . which includes 
the following theoretical codes. For the RGE 
running of the soft SUSY-breaking parameters, 
it uses SoftSUSY [14], which is combined con- 
sistently with the codes used for the various 
low-energy observables: FeynHiggs [15HT8] is 
used for the evaluation of the Higgs masses and 
a u USY (see also [T9ll20] ). for the other electroweak 
precision data we have included a code based 
on [HE], SuFla [2l][22] and Superlso [23 | f24 ] are 
used for flavor-related observables, and for dark- 
matter-related observables MicrOMEGAs [33] and 
DarkSUSY [2S1 are included. In the combination 



of the various codes, MasterCode makes extensive 
use of the SUSY Les Houches Accord [27 1 128 ] , 

The experimental constraints used in our anal- 
yses arc listed in Table 1 of Ref. [T]. One im- 
portant comment concerns our implementation of 
the LEP constraint on Mh- The value quoted in 
the Table 1 of Ref. PQ, M H > 114.4 GcV, was 
derived within the SM [29], and is applicable to 
the CMSSM, in which the relevant Higgs cou- 
plings are very similar to those in the SM [3U1I3TJ , 
so that the SM exclusion results can be used, 
supplemented with an additional theoretical un- 
certainty: we evaluate the x 2 {Mh) contribution 
within the CMSSM using the formula 
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with M h imit 



115.0 GeV for M h < 115.0 GeV. 
Larger masses do not receive a x {Mh) contribu- 
tion. We use 115.0 GeV so as to take into account 
experimental systematic effects. The 1.5 GeV in 
the denominator corresponds to a convolution of 
the likelihood function with a Gaussian function, 
<i>i.5(a:), normalized to unity and centered around 
Mh, whose width is 1.5 GeV, representing the 
theory uncertainty on Mh [17j . In this way, a 
theoretical uncertainty of up to 3 GeV is assigned 
for ~ 95% of all Mh values corresponding to one 
CMSSM parameter point. The 1.1 GeV term in 
the denominator corresponds to a parameteriza- 
tion of the CL S curve given in the final SM LEP 
Higgs result [2S|. Within the NUHM1 the situ- 
ation is somewhat more involved, since, for in- 
stance, a strong suppression of the ZZh coupling 
can occur, invalidating the SM exclusion bounds. 
In order to find a more reliable 95% C.L. exclu- 
sion limit for Mh in the case that the SM limit 
cannot be applied, we use a specially adopted pro- 
cedure outlined in Ref. [T] (using the results of 
the code HiggsBounds [32] that incorporates the 
LEP [33] (and Tevatron) limits on neutral Higgs 
boson searches). 

Our statistical treatment of the CMSSM and 
NUHM1 makes use of a large sample of points 
(about 2.5 x 10 7 ) in the SUSY parameter spaces 
obtained with the Markov Chain Monte Carlo 
(MCMC) technique. Our analysis is entirely fre- 
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Figure 1. The (mo,m 1/2 ) plane in the CMSSM (left) and the NUHM1 (right). The dark shaded area at 
low mo and high TO1/2 is excluded due to a scalar tau LSP, the light shaded areas at low TO1/2 do not 
exhibit electroweak symmetry breaking. Shown in both plots are the best-fit point, indicated by a filled 
circle, and the 68 (95)% C.L. contours from our fit as dark grey /blue (light grey /red) overlays, scanned 
over all tan/3 and A$ values. The 95% C.L. exclusion curves (hadronic search channel) at CMS with 
1 (0.1) fb _1 at 7 TeV center-of-mass energy is shown as green/light gray (black) solid curve. 



qucntist, and avoids any ambiguity associated 
with the choices of Bayesian priors. 



3. PREDICTIONS 

MASSES 



FOR SPARTICLE 



For the parameters of the best-fit CMSSM 
point wc find m = 60 GcV, mi/2 = 310 GeV, 
A = 130 GcV, tan/3 = 11 and \x = 400 GcV, 
yielding the overall x 2 /N dof = 20.6/19 (36% 
probability) and nominally M} % = 114.2 GeV. 
The corresponding parameters of the best-fit 
NUHM1 point are m = 150 GeV, m 1/2 = 
270 GeV, A = -1300 GeV, tan/3 = 11 and 
m\ = m 2 h2 = —1.2 x 10 6 GeV 2 or, equivalently, 
H = 1140 GeV, yielding x 2 = 18-4 (correspond- 
ing to a similar fit probability as in the CMSSM) 
and M h = 120.7 GeV. 

In Fig. Q] we display the best- fit value and the 
68% and 95% likelihood contours for the CMSSM 
(left plot) and the NUHM1 (right plot) in the 
{mo, mi/ 2 ) plane, obtained as described above 
from a fit taking into account all experimental 
constraints. We also show exclusion contours for 
the hadronic search mode (jets plus missing en- 
ergy) at CMS. The green (light gray) solid line 



shows the 95% C.L. exclusion contour for CMS for 
1 ftT 1 at y/s = 7 TeV [34]. The black solid line 
shows the corresponding results for only 0.1 fLV . 
(Similar results hold for ATLAS.) One can see 
that with 1 ftV 1 the best-fit points can be tested, 
together with a sizable part of the whole 68% C.L. 
preferred regions. In the case of the NUHM1 
(right plot) nearly the 68% C.L. region could be 
probed. 

Wc now review the results for the predic- 
tions of sparticles masses in the CMSSM and the 
NUHM1, which are summarized in Fig.[2j The re- 
sults for the CMSSM spectrum are shown in the 
left plot, and for the NUHM1 in the right plot. 
We start our discussion with the gluino mass, m§. 
In both the CMSSM and the NUHM1, the best-fit 
points have relatively low values of rrig ~ 750 and 
~ 600 GeV, respectively. These favored values 
are well within the range even of the early opera- 
tions of the LHC with reduced centre-of-mass en- 
ergy and limited luminosity, in accordance with 
the result of Fig. [TJ However, even quite large 
values of irig < 2.5 TeV are allowed at the 3- 
er (Ax 2 = 9) level (not shown in Fig. HJ. The 
LHC should be able to discover a gluino with 
m g ~ 2.5 TeV with 100/fb of integrated lumi- 
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Figure 2. Spectra in the CMSSM (left) and the NUHM1 (right) [J. The vertical solid lines indicate the 
best-fit values, the horizontal solid lines are the 68% C.L. ranges, and the horizontal dashed lines are the 
95% C.L. ranges for the indicated mass parameters. 



nosity at -y/s = 14 TcV ;35,|3E], and the pro- 
posed SLHC luminosity upgrade to 1000/fb of in- 
tegrated luminosity at ^/s = 14 TeV should per- 
mit the discovery of a gluino with nig ~ 3 TeV 
[37]. 

The central values of the masses of the su- 
persymmetric partners of the u, d, s, c, b quarks 
are slightly lighter than the gluino, as seen in 
Fig. [2] The difference between the gluino and the 
squark masses is sensitive primarily to m$. The 
SUSY partners of the left-handed components of 
the four lightest quarks, the gx, are predicted to 
be slightly heavier than the corresponding right- 
handed squarks, qn, as seen by comparing the 
mass ranges in Fig. [2] As in the case of the gluino, 
squark masses up to ~ 2.5 TeV are allowed at the 
3-cr level. Comparing the left and right panels, we 



sec that the squarks are predicted to be somewhat 
lighter in the NUHM1 than in the CMSSM. 

The scalar taus as well as the other scalar lep- 
tons are expected to be relatively light, as can be 
seen in Fig. [2j They would partially be in the 
reach of the ILC(500) (i.e. with y/s = 500 GeV) 
and at the 95% C.L. nearly all be in the reach 
of the ILC(IOOO) [38j[39]. This also holds for the 
two lighter neutralinos and the light chargino. 

4. PREDICTIONS FOR THE LIGHTEST 
HIGGS MASS 

Finally we discuss the likelihood functions for 
M h within the CMSSM and NUHM1 frameworks 
obtained when dropping the contribution to \ 2 
from the direct Higgs searches at LEP. The results 
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Figure 3. The \ 2 functions for M h in the CMSSM (left) and the NUHM1 (right) [I], including the 
theoretical uncertainties (red bands). Overlayed is the blue band curve for the SM Higgs boson taken 
from Ref. [40]. Also shown is the mass range excluded for a SM-like Higgs boson (yellow shading). The 
LEP limits are not included in the fit. 



are shown in the left and right panels of Fig. [3J 
respectively, together with the corresponding SM 
result g0]. 

It is well known that the central value of the 
Higgs mass in a SM fit to the precision elec- 
troweak data lies below 100 GeV [40], but the 
theoretical (blue bands in Fig. [3J and experimen- 
tal uncertainties in the SM fit are such that they 
are still compatible at the 95% C.L. with the di- 
rect lower limit of 114.4 GeV [29] derived from 
searches at LEP. In the case of the CMSSM and 
NUHM1, one may predict Mh on the basis of the 
underlying model parameters, with a l-o~ uncer- 
tainty of 1.5 GeV [17], shown as a red band in 
Fig. [3j Also shown in Fig. [3J are the LEP exclu- 
sion on a SM Higgs [29] and the Tevatron exclu- 
sion around Mh ~ 165 GeV glj (yellow shad- 
ing). The LEP exclusion is directly applicable to 
the CMSSM, since the h couplings arc essentially 
indistinguishable from those of the SM Higgs bo- 
son [30l[3T] , but this is not necessarily the case in 
the NUHM1, as discussed earlier. 

In the case of the CMSSM, we see in the left 
panel of Fig. [3J that the minimum of the x 2 func- 
tion occurs below the LEP exclusion limit. The 
fit result is still compatible at the 95% C.L. with 
the search limit, similarly to the SM case. As 
discussed above, a global fit including the LEP 



constraint has acceptable x 2 - In the case of the 
NUHM1, shown in the right panel of Fig. [3J we 
see that the minimum of the x 2 function oc- 
curs above the LEP lower limit on the mass of 
a SM Higgs at M h m 121 GeV. Thus, within the 
NUHM1 the combination of all other experimen- 
tal constraints naturally evades the LEP Higgs 
constraints, and no tension between Mh and the 
experimental bounds exists. 
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